Ground State Magnetization of Polymerized Spin Chains 



Daniel C. Cabra 1,2 and Marcelo D. Grynberg 1 
1 Departamento de Fisica, Universidad Nacional de la Plata, C.C. 67, (1900) La Plata, Argentina. 
Facultad de Ingenieria, Universidad Nacional de Lomas de Zamora, 
Cno. de Cintura y Juan XXIII, (1832), Lomas de Zamora, Argentina 

We investigate the ground state magnetization plateaus appearing in spin 1 /2 polymerized Heisen- 
berg chains under external magnetic fields. The associated fractional quantization scenario and the 
exponents which characterize the opening of gapful excitations are analyzed by means of abelian 
bosonization methods. Our conclusions are fully supported by the extrapolated results obtained 
from Lanczos diagonalizations of finite systems. 
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Models of low-dimensional magnets, such as strongly 
correlated quantum spin chains yj] and ladders are 
currently receiving renewed and systematic attention for 
a variety of reasons. Amongst the most remarkable are 
the spin-Peierls dimerization instability ||,0], the Hal- 
dane conjecture j^], and novel concepts such as fractional 
quantization and topological energy gaps (l]||. These 
rather complex phenomena are largely owing to quan- 
tum fluctuations of individual spins which tend to restore 
the rotational symmetry of the ground state. Depending 
on the exchange interactions, fluctuations can manifest 
themselves collectively into many possible ground states, 
particularly in lower dimensions where their effects are 
most severe 0. 

A wealth of issues have been addressed experimentally 
to confirm these expectations in a series of quasi-one- 
dimensional compounds ||. After a vast body of re- 
search, it is by now well established that half-integer- 
spin chains are massless whereas integer ones are gapful 
(see also ||). In spite of the availability of a number of 
excellent realizations of one dimensional Heisenberg anti- 
ferromagnets |^,^| , detailed measurements of spin excita- 
tions however, have remained confined within the limits 
of applied magnetic fields which are low with respect to 
the exchange interactions thus, precluding a comparison 
with theory. Nevertheless, advanced neutron scattering 
experiments spanning large magnetic fields regions in rel- 
atively low exchange coupling materials, e.g. CS2CUCI4, 
are now beginning to appear (see e.g. |10(| ). 

Magnetization properties of spin chains have been the 
subject of intensive investigations for quite a while. Im- 
portant work includes that of Parkinson and Bonner p| , 
and has recently gained new interest due to the appear- 
ance of magnetization plateaus studied by several authors 
in a variety of systems [^2|- Jig] . As a contribution to the 
understanding of massive spin excitations in high mag- 
netic fields, here we consider an alternative scenario for 
their low-temperature realization. Specifically, we study 
the interplay between explicit breaking of full transla- 
tional symmetry and applied magnetic fields, say along 
the z-direction, in spin-1/2 Heisenberg chains at T = 0. 
This can be conveniently described by a set of Heisenberg 
antiferromagnets in which the exchange coupling interac- 



tions J n are all equal but one every p sites, i.e. 



J(i 

J, 



n/pe Z , 
otherwise . 
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It should be stressed that periodic arrays of couplings 
are relevant to the study of ferrimagnetic materials pToj ] 
and that also one dimensional dimcrized and trimerized 
materials are known to exist PUl2|]. 

The Hamiltonians of our polymerized chains in the 
presence of a (dimensionless) magnetic field h applied 
along the z-axis are thus given by: 



n=l 



■ ^2 S " 

71=1 
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Here, the S n are spin-1/2 operators, whereas periodic 
boundary conditions are assumed along the L sites of 
the chain £ Z). Despite their simplicity, we will 
show however that these Hamiltonians entail a highly 
non-trivial magnetic behavior controlled solely by the 
chain periodicity p a nd the external field h. 

As is known |20[ . the full translationally invariant 
(FTI) S = h chain remains gapless for all magnetizations 
(M) = -| (Sn^n)) U P to a saturation field where each 
of the L individual spins becomes fully polarized. On 
general grounds however, the Lieb-Schultz-Mattis theo- 
rem |l4||2~l| ] indicates that FTI Hamiltonians of arbitrary 
spin-5, can be gapful provided the magnetization per 
spin (M) satisfies (5* — \{M)) e Z. Such gapful exci- 
tations should be reflected through the presence of mag- 
netization plateaus, in principle at these special values 
of (M). However, notice that the above theorem does 
not prove the existence of this quantization scenario as it 
refers to non-magnetic excitations, i.e. modes which pre- 
serve the total magnetization. Nevertheless, magnetiza- 
tion plateaus have been extensively conjectured and ob- 
served in both Peierls dimerized (|,[l^] an d trimerized JlJ] 
spin chains, as well as in frustrated p/flJlSfl , anisotropic 
spin- 1 systems 1 13 and ladders models~i|16(|. They all are 
examples of a rather subtle phenomenon namely, frac- 
tional quantization of a macroscopic physical quantity 
under external varying fields. Here, we examine this sit- 
uation for the whole class of non- homogeneous chains M), 



1 



attempting to extend and systematize aspects of quanti- 
zation emphasized in those studies. 

Following a recent analysis discussed as in ]l5| , ^6| , ^8| 
and comprehensively accounted in PJ^, we will apply 
the by now standard method of abelian bosonization to 
Eq. (g). It is well known that the low-energy properties 
of the FTI Hcisenberg chain, (8 = 0), are described by a 
c = 1 conformal field theory of a free bosonic field com- 
pactified at radius R for any given magnetization (M) 
(see e.g. |23] ). The functional dependence of R can be 
obtained from the exact Bethe Ansatz solution by solv- 
ing a set of differential equations obtained in [|3and |2j| 
(for a fuller derivation consult for instance Ref. Jl6[ ). Ex- 
ploiting this knowledge, the bosonized expression of the 
low-energy effective Hamiltonian (||) in the homogeneous 
case (5 = 0, can be readily shown to adopt the form 



H 



dx 



^{u\x)+R\{M)){d x <P{x)f} , (3) 



with II = ^d x 4>, and 4> = 4>l+4>r, 4> = 4>L-<t>R- Here, the 
effect of the magnetic field h enters through the radius of 
compactification R((M)). This radius governs the con- 
formal dimensions, in particular the conformal dimension 

of a vertex operator e ll3< ^ is given by ^^J^ ■ Within the 

framework of the theory of Luttinger liquids, it is worth 
pointing out also that the compactification radius is re- 
lated to the parameter K by R 2 = J^. 

The bosonized expressions for the spin operators read: 



S!(x) 



and 



2tt dx 



+ const. : cos(2kpX + v47T0i) : + 



(1 + const. : cos(2kpX + V47T0 



(M) 
2 

(4) 
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where the colons denote normal ordering with respect 
to the ground state with magnetization (M). Now we 
apply this methodology to compute the effective form of 
the interaction. After some algebra, in the limit of weak 
polymerization 6 <C 1, it can be readily shown that the 
most relevant perturbation term is given by 



L/p 



H int » 8 ^2 cos(2k F (px' + h + V^ttc 



(6) 



This operator will survive in passing from the lattice 
to the continuum model, assuming that the fields vary 
slowly, only when the oscillating factor exp(i2px'kp) 
equals one. Since the Fermi level is given by kp = 
■|(1 — (A/)), this in turn will happen when the condition 



l(i-(M))eZ 



(7) 



is satisfied . 

We can now study when a plateau will appear in the 
magnetization curve for the polymerized chain (Q). To 
do that, we first have to see which are the values for 
the magnetization where there could be a plateau, for a 
given value of the period p, (i.e. solve for (0)) and then 
we need to evaluate the scaling dimension of the operator 
(||), which at zero loop is given by 



d = 



1 



AttR 2 '' 



(8) 



which is in turn governed by the radius of compactifi- 
cation as we already stressed. By virtue of the lower 
bound of the compactification radius |24|,|25|], namely 
R((M)) > R(±l) = 1/(2 VtF) , h follows from Eq. (g) 
that d < 2 for all magnetizations \(M)\ < 1. This en- 
sures the relevant character of the operator (g) which in 
turn survives in the continuum limit whenever Eq. (j?]) is 
satisfied. Therefore, we can conclude that constraint (Q) 
constitutes a sufficient condition ultimately responsible 
for the appearance of magnetization plateaus and mas- 
sive spin excitations. This is our main result. 

We now turn to a numerical finite-size analysis. In 
Figs. l(a)-(e) we display a variety of magnetization 
regimes as a function of both polymerization parameters 
J' /J = 1 — 5 and applied magnetic fields h. This is a 
rather compact form of representing conventional mag- 
netization curves for different polymerization strengths. 
Here, each line is associated to successive values of (M) 
which decrease monotonically from top to bottom, as 
they should for a non- frustrated system. The results 
were obtained from exact diagonalization of finite sys- 
tems via a recursion type Lanczos algorithm [ p6| applied 
on each magnetization subspace S z = {0, 1, ... ,L/2}. To 
avoid the formation of spurious interfaces, even multiple 
lengths of the lattice periodicity were taken throughout. 
Using fully isotropic chains up to L = 24 sites with pe- 
riodic boundary conditions, our numerical analysis sup- 
ports entirely the quantization constraint (0). 

As expected, the ground state 'phase diagrams' ex- 
hibit bands of empty regions corresponding to the mag- 
netization plateaus of S z referred to above, while regions 
filled with magnetization lines reflect smooth magnon ex- 
citations arising in the thermodynamic limit L — > 00. 
It can be readily observed that for chains of periodic- 
ity p > 1 (dimers, trimers, etc), a plateau- like struc- 
ture emerges precisely at the rational magnon densities 
(M) = 1 — 2q/p, (q — 0, 1, ...p), implicit in the general 
scenario of Eq. (J?]) . It is worth remarking on the robust- 
ness of this topological constraint as similar results con- 
tinue to hold for anisotropic (XXZ) chains, the plateaus 
always appearing at the same values of (M ) . 

We can also predict the behavior of the mass gap 
(width of the plateau) , with the polymerization strength 
8 by means of a simple zero loop computation J27]] . Aside 
logarithmic corrections to the case (M) = 0, this yields 



(9) 
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with d given as in (g). 

To enable an independent check of this result, we now 
turn to the issue of extrapolating the numerical finite size 
estimates of the mass gaps gL towards their correspond- 
ing thermodynamic limits. Note, on one hand, that any 
extrapolation procedure by necessity assumes that the 
asymptotic behavior applies to the values of L within 
reach. However, it is known [ pp[ that finite size correc- 
tions to the gap in the excitation spectrum of the homo- 
geneous Heisenberg chain vary slowly as In (lni)/ln 2 (i), 
thus affecting the results over a wide range of sizes. In 
fact, as can be seen in Figs. (||), this turns out to be the 
case also for weak polymerization regimes, 8 — > 0, where 
finite size effects are more pronounced. Therefore, in 
studying numerically the mass gap behavior obtained in 
Eq. (M) we are confronted to restrict considerations to the 
non-critical region < \S\ < 1, yet suitable to test inde- 
pendently the correctness of our bosonization approach. 

To estimate the actual masses in the limit L — > oo, 
we fitted the whole set of finite-size results (even inte- 
ger multiples of p within the range 4 < L < 24), using 
both linear, and logarithmic type methodologies of con- 
vergence pq] , i.e. 



gL ^ g 



Ae 



-BL 



g L ~ g + A/L 1 



(10) 



Either extrapolation procedure yields basically the same 
result with at least 3 significant digits. This latter vari- 
ation ultimately gives an estimative idea of the lower 
bound of the extrapolation error. The reliability of our 
results was checked also by comparing the trend arising 
from smaller systems (L < 20). When the critical re- 
gion is approached however, the accuracy differs widely, 
particularly for \S\ < 0.2. 

Although there are alternative extrapolation algo- 
rithms which do not involve fits to specific forms |2q | 
we should hasten to add however, that their efficiency 
depends strongly on the abundance of data. In our case, 
this is translated in the availability of matching sizes, 
already constrained by both the periodicity p and the 
antifcrromagnetism. Nevertheless, we were able to find a 
remarkable agreement with the compactification radius 
comprehended in Eq. (@) and the exponents of Eq. (||). 
The results are shown in Fig. 2 where we display respec- 
tively the gap openings around (M) =0, 1/3, 1/2, for 
p = 2, 3 , and 4. The dimerized case reproduces the well 
known 2/3 exponent predicted in || and corroborated 
subsequently by diverse numerical studies 29 1. To our 



knowledge however, opening exponents forp > 3 (see Fig. 
2), have not been elucidated yet by other investigations. 

Finally, it is instructive to comment further on the role 
of quantum fluctuations namely, the tendency of spins to 
spontaneously tilt occasionally due to the Heisenberg un- 
certainty relations, and their relevance to our results 0. 
For classical spins, e.g. Ising and n-vector models, the 
interplay between dimensionality and statistical fluctua- 
tions, though crucial in determining phase transitions, is 
not sufficient to entail the fractional behavior studied so 



far. In fact, an elementary transfer matrix calculation 
shows that the Ising equivalent of (0) wipes out all but 
two magnetization plateaus, namely (M) = 0, 2/p, (even 
p > 2), their widths behaving linearly with 5. Thus, it 
is worth pointing out that Eqs. (Q) and (Q), in contrast, 
constitute a genuine macroscopic quantum effect. 

In summary, we have presented a bosonization picture 
that accounts for the fractional quantization observed in 
a class of non-homogeneous Heisenberg antiferromagnets. 
All low-energy exponents which characterize the opening 
of gapful excitations have been obtained and treated on 
an equal footing while checked with Lanczos diagonaliza- 
tions of finite systems. Aside these theoretical pursuits, 
we trust our study will help to convey a clearer under- 
standing of the many characteristics present in real low 
dimensional magnets. A similar analysis in polymerized 
ladder systems is in progress. 
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FIG. 1. Magnetization contours of finite polymerized 
chains for (a) p = 2, L = 24, 20, 16 ; (b) p = 3, L = 24, 18, 12 ; 
(c) p = 4, L = 24, 20, 16 ; (d) p = 5, L = 20, 10 (full and 
dotted lines respectively) and; (e) p = 6, L = 24, 18, 12 . Ex- 
cept for (d), full, dashed, and dotted lines stand respectively 
for large, medium and small sizes. They denote all accessible 
magnetizations, whereas their values decrease from top to bot- 
tom. Though numerical accuracy in h/ J is bounded by 10~ 7 , 
size effects become evident for J' = J, as no plateaus (empty 
wide bands), should occur in the thermodynamic limit. 
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FIG. 2. Extrapolated values of the gap for p — 2 around 
(M) = 0, p = 3 for (M) = 1/3 and, p = 4 with (M) = 0, 1/2. 
Solid lines are guide to the eye whereas slopes of dashed lines 
denote the estimated opening exponents, namely (in descend- 
ing order), 0.8(1), 0.77(10) and, 0.66(10), ((M) = 0). To 
improve the clarity of the figure, the uppermost curve was 
shifted multiplying the gap by a scale factor 2.5. 
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